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Abstract 

The supersymmetric technique is applied to computing the average spectral density near zero energy 
in the large-A'^ limit of the random-matrix ensembles with zero eigenvalues: B, Dlll-odd, and the chiral 
ensembles (classes j4III, BD\, and CII). The supersymmetric calculations reproduce the existing results 
obtained by other methods. The effect of zero eigenvalues may be interpreted as reducing the symmetry of 
the zero-energy supersymmetric action by breaking a certain abelian symmetry. 

There exists a remarkable correspondence between large families of random-matrix ensembles and symmetric 
superspaces. It has been shown by Zirnbauer that in the large- limit [N is the matrix dimension) correlation 
functions in random-matrix ensembles may be represented as integrals over appropriate Riemannian symmetric 
superspaces (with dimensions independent of N) [1]. This relation to symmetric superspaces is based on 
Efetov's supersymmetric technique introducing auxiliary anticommuting (Grassmann) variables in order to 
directly average correlation functions over the statistical ensemble [2]. 

At the same time, the random-matrix ensembles are known to be in one-to-one correspondence with sym- 
metric spaces (Cartan symmetry classes) [3, 4]. The classification of Zirnbauer thus establishes a correspondence 
between large families of symmetric spaces and Riemannian symmetric superspaces [1]. The random-matrix en- 
sembles with zero eigenvalues were not included in the original classification, and later it became apparent that 
the zero eigenvalues in random-matrix ensembles are related to the reducibility of the correseponding symmetric 
superspaces [5, 6]. 

In this paper, I study this relation by explicitly calculating the average density of states in all random-matrix 
ensembles with zero eigenvalues. There are five such ensembles (Table 1): class B [so(A^) matrices at odd iV], 
class DIITodd [so(2A^)/u(A^) matrices at odd N], and the three chiral ensembles: unitary AIII, orthogonal BDl, 
and symplectic CII. In a physical context, the ensembles B and Z?III-odd appear in vortices in superconductors 
with odd pairing [5, 7], the chiral classes — in QCD [8, 9]. The zero levels in these ensembles occur as a 



Table 1: Random- matrix ensembles with zero eigenvalues 



Cartan 
class 


Symmetric space 
(compact type) 


/3 


Number of zero 
a . , 

eigenvalues m 


B-D 


SO(iV) 


2 


2m 


m=0, even N 


Dill 


S0{2N)/13{N) 


4 


l-|-4m 


m—l, odd N 


Alll 


SU(p-}-q)/S(U(p)xU(g)) 


2 


l+2m ' 




BDl 


SO(p+g)/SO(p)xSO(q) 


1 


ra 


> m=\p-q\ 


CII 


Sp(p+g)/Sp(p)xSp(q) 


4 


3+4m 
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consequence of the symmetry inverting energy (E —E) combined with the odd dimension (for classes B and 
Dlll-odd) or with the dimensional imbalance between the two chiral sectors (for the chiral classes). Table 1 
also lists the values of the parameters a and (3 appearing in the joint probability distribution for energy levels 

LOi'. 

dP{uJi, . . . , a;„) oc n - ^jf n '^i'^i (1) 

i<j i 

{(3 determines the strength of repulsion between levels, a — the strength of repulsion from zero). 

Previously, a supersymmetric calculation of the misroscopic spectral density for the chiral unitary case was 
done in Ref. [6], and the case of class B was studied in Ref. [5] (in the context of class BD, which is the average 
of classes B and D) . I include these cases for completeness in the corresponding sections. 

Prom the supersymmetric calculations for the five random-matrix ensembles, we find that zero levels in 
random-matrix ensembles manifest themselves in reducing the symmetry of the supersymmetric action at zero 
energy. In the absence of zero levels, this action (a function of the supermatrix Q in Efetov's technique [2]) is 
invariant with respect to the full supergroup preserving the linear constraints on Q (latter being determined 
by the symmetries of random matrices). For ensembles with zero levels, the zero-energy action is invariant 
with respect to only a normal subgroup of this supergroup, but breaks the remaining abelian symmetry. In the 
large- A?^ limit, the integral over Q is dominated by the saddle-point manifold. This manifold is a Riemannian 
symmetric superspace [1], and for ensembles admitting zero levels it is not irreducible: it may be split into 
orbits of the normal subgroup of the full symmetry (super) group. The quotient by this normal subgroup is an 
abelian (conventional, not super) group (Z2 for classes B-D and Dili, and GL(1) for the chiral classes). If the 
random-matrix ensemble contains zero levels, the action is not invariant with respect to this residual abelian 
group, but transforms according to one of its one-dimensional representations. 

The paper is organized as follows. In the next section, I review the results for the average spectral density in 
the random-matrix ensembles with zero levels. Next, I describe the details of the supersymmetric calculations 
for each of the five random-matrix ensembles. The calculation for the ensemble B-D is presented in somewhat 
more detail, and in the subsequent sections the repeating steps of the derivations are described only briefly. In 
the last section I discuss common features of these calculations specific for ensembles with zero levels. 



1 Spectral density in random-matrix ensembles with zero levels 

In this section I review the results for the average spectral density in the vicinity of the zero eigenvalue. All 
these results are known and were previously derived by other methods. The reader may use this section as a 
quick reference. 

In what follows we consider zero-curvature random-matrix ensembles and treat them as quantum-mechanical 
Hamiltonians. Accordingly we use quantum- mechanical terminology such as "energy levels", "inter- level spac- 
ing", etc. 

In an ensemble of random matrices of size N , with a fixed dispersion of matrix elements, the inter-level 
spacing in the middle of the spectrum scales as N'^/"^ for large A''. If we measure the energy E in the units of 
this inter-level spacing A, the correlation functions in the vicinity of zero energy (middle of the spectrum) have 
a finite and universal limit as A'' — > 00. In this paper I am interested in the average density of states p{x) as 
a function of dimensionless energy x = E/A. This function gives the average number of energy levels in any 
interval [a;b]: 

(»^>[a;b] / p{x)dx (2) 

(In the case of symplcctic ensembles _DIII and CII, all energy levels are doubly degenerate, and for counting 
purposes every degenerate pair of states will be counted as a single level). The function p{x) is symmetric 
p{x) = p{—x) and has the normalization \\-mx^oo p(x) = 1. For an ensemble with m zero levels, p{x) has a 
(5-functional contribution at x = 0: p{x) = mS{x) + p{x), where p{x) is continuous at a; = 0. 

The results for the average density of states p{x) in the ensembles studied in this paper are the following 
(defining y = 2Tr\x\, m is the number of zero levels in the chiral ensembles). 
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Class D: 



Class B: 



1 + 



siny 



Class £)III-even: 

Class £>III-odd: 

Class All! (chiral unitary): 



siny ^, . 
1 + 5{x) 

y 



[J[iy)My) + JKy)]+^Ji{y) 



|y [J[iy)My) + JUy)] - ^Mv) + Hx) 



Jmi^) - Jm-l{^)Jm+l{^) + m6{x) 



Class BDl (chiral orthogonal): 

f (f ['^-(f)-'^— i(f)'^'»+i(|)] +-^'"(|)^'"(|)) +"^^(^) 
Class CII (chiral symplectic): 

I [y [Jimiy) - J2m-l{y)J2m+l{y)\ " J2m{y)R2m{y)) + m5{x) 

where the functions i?„ and i?„ are defined as: 

R 



:„{Z) = 1 - Rn{z) = [ Jniz')dz'. 

Jo 



(3) 
(4) 
(5) 
(6) 

(7) 
(8) 
(9) 
(10) 



These results were previously derived by other methods. The results (3) and (4) are presented in the book 
of Mehta [10]. They are also straightforward to obtain from mapping of level statistics onto free fermions. A 
supcrsymmctric approach to classes B and D was developed in Rcf. [5]. The result (5) was found by Nagao 
and Slevin [11] and by Altland and Zirnbauer [12] (contrary to their claim, their result is identical to the result 
of Nagao and Slevin after some algebraic manipulations with Bessel functions). The result (7) was obtained in 
the works of Verbaarschot and Zahcd [13], Nagao and Slevin [14], and Forrester [15]. Also, a supcrsymmctric 
calculation of (7) at m = was reported in [16], and then at arbitrary m in [6]. To make this paper self- 
contained, I repeat the derivation of Ref. [6] in the corresponding section. The particular case of the formula 
(8) at m = 1 can be foimd in [11]. The case of arbitrary m was treated in [17] and [18]. The latter work also 
contains the answer for the ensemble CII. The results of [17] and [18] are presented in the form of complicated 
integrals. The simple formulas (8) and (9) were later reported in Refs. [19, 20, 21]. 

The average spectral densities (3)-(9) are plotted in Fig. 1. 

2 Remarks about notation, supergroups and superspaces 

In this section I explain some notational conventions used in the subsequent sections. The calculations involve 



supermatrices acting in a superspace which has the structure C ig) C or 



depending on the 



symmetry class. One of the factors refers to the Fermi-Bose (FB) sectors and defines the supersymmetric 

grading. The one or two remaining factors arc cither produced by additionally doubling the dimension to take 
into account the symmetries of the random- matrix ensemble (in classes B-D, Dill, and BDl) or are originally 
present in the matrix structure in the random-matrix ensemble (in class Dili and in the three chiral classes). 
These will be labeled to as "particlc-holc" (PH) or "1-2" sectors, without stressing the physical meaning of 
this terminology. In supermatrices, the FB sectors will be graphically divided by solid lines [see, for example. 
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eq. (12)], with the BB sector in the upper left corner, and the FF sector in the lower right corner. When the 
matrices also act in the PH or 1-2 spaces, in order to avoid confusion, these spaces will be explicitly mentioned 
as a subscript, from the outermost division to the innermost subdivision [see, for example, eqs. (58)-(60)]. The 
empty spaces in matrices denote zeroes. 

The supergroups GL(n|m) and 0S'p{n\'2rn) appearing in our supersymmetric constructions are defined as 
follows. The complex supergroup GL(n|m) consists of all invertible supermatrices of dimension n + m. The 
complex supergroup 0Sp(n|2m) is the subgroup of GL(n|2m) obeying the relation 

5"' = 75V\ (11) 

where 



/ In 






Im 




-Iru 



(12) 



Its support is the direct product of 0(n) and Sp(m). More details about these supergroups may be found in 
Refs. [22, 1]. The reader may also refer to Ref. [2] for conventions regarding manipulations with supermatrices. 
The Lie superalgebras of 0Sp(n|2m) and GL(n|m) are denoted as osp(n|2m) and gl(n|m). 

To distinguish between fermionic and bosonic sectors [which is important when performing integration, 
either compact or non-compact, see below], we shall reserve the notation 0Sp(nl2m) for the supergroup with 
0(/i) in the bosonic and Sp(m) in the fermionic sector. The same supergroup with 0(n) in the fermionic and 
Sp(m) in the bosonic sector we denote as SpO(n]2m). Also, we use the notation SpS0(n]2m) for the connected 
component of SpO(nl2m) [with the unit superdeterminant] . 

The notation Sp(to) in this paper refers to the symplectic group of 2to x 2m matrices. This notation agrees 
with Refs. [1, 3, 7], but differs from Ref. [22] where the same group is denoted Sp(2m). 



3 Classes B and D 

In this section we use the supersymmetric technique to compute the density of states for the so(A^) random 
matrices (class D for even TV, class B for odd A^). This is the simplest of the five examples considered in this 
paper, and we describe it in more detail to demonstrate the technique of the calculation. The calculation follows 
the prescription described in detail by Zirnbauer [1]. 

The random- matrix ensembles B-D is unitary (with (3 = 2). The supermatrix Q used in the calculation 
of the average density of states has dimension 2 + 2 (2 bosonic and 2 fermionic dimensions), is parameterized 
by 4 -h 4 independent variables and is an element of osp(2]2) Lie superalgebra. The saddle-point manifold 
has dimension 2 -|- 2, and thus the density of states in the large- A limit is computed as an integral over two 
commuting and two Grassmann variables. 

The random- matrix ensemble B-D consists of purely imaginary antisymmetric matrices H: 

Hab = = -Hba, a,b=l,...,N. (13) 
The matrix elements have independent Gaussian distributions: 

dP{H) o, H exp (-^1^) dHab, (14) 

a>b ^ ' 

SO that the averages of any number of matrix elements arc given by the Wick rule together with the pair average 

{HabHa'b') = V^{5ab'Sba' " Saa'^bb')- (15) 

Prom the calculation below we shall see that the energy unit defined as 

A = ^ (16) 
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plays the role of the average level spacing near zero energy. 

The average density of states can be found by differentiating the generating function 



J dP{H) 



det{H - wpA) 
det {H -ojb^Y 



(17) 



where the integration is performed over the ensemble of random matrices H] Ub and wp are auxiliary variables 
(complex numbers). We included the energy scale A in the definition (17) to make ujb and ujp dimensionless. 

The ensembles considered in this paper have an E ^ —E symmetry, which leads to the symmetry of the 
generating function: 

Z{ujb,ujf) = Z{-u;b,-uJf) = (-1)'"^(-wb, c^f), (18) 

where to is the number of zero levels. In the supersymmctric calculation below we neglect the overall sign of 
Z{u!b,ujf), but restore it at the end of the calculation from the condition Z{lu,lu) = 1 and from positiveness of 
the density of states. 

The two determinants in (17) may be written as Gaussian integrals over bosonic and fermionic variables 
(auxihary fields). Introducing the (A^ + 7V)-component supcrvcctor tpa — {ipBa,fpFa),o, = 1,...,N, and per- 
forming the integration over dP{H), we arrive at the partition function for interacting superfields [the common 
energy scale v drops out already at this step, thanks to our including A in (17)]: 



Z{ujb,ujf) 



£»(V'^^/')exp (-^^Via^/xa - I [(V'iaV'Mfe)(^!bV'i'a) - {i^iai^tib){i^laM 



(19) 



where /i, i' are fcrmion-boson indices. In the integral, the Grassmann components in ip and ■0^ are treated as 
independent variables (total 2N Grassmann variables). The integral over bosonic components of ip and tp^ is 
taken over the 2A''-dimensional real submanifold (V'^)sa = {''pBa)*- 

To decouple the interaction with the Q-matrix, it is necessary to double the dimension of vector i/j. Combine 
the old superfields tp^a into the new ones: 




SO that 
where 



^ n. — 



/ 1 

-1 



V 



-IpB 



1 

1 0/ 



(20) 



(21) 



(22) 



In terms of the supervectors ^ and the partition function may be rewritten as 

1 



Z{uJb,U1f 



£)* exp 1 -2STr 



where 



ujb 






U)F 



1 
-1 



(23) 



(24) 



(With this definition of ui we in fact change the sign oiup, which may result in the change of sign of Z{u)b, i^f), 
according to (18). We shall not control the overall sign of Z{uJB,ijJF), but restore the correct sign at the end of 
the calculation.) 
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The matrix ("fa^a) has the exphcit form 



/ Qb 
X 



-X 
-Qb 



\ 



a 
P 



-P 
-a 



a 
P 



P 
a 



\ 



Qf 






-Qf j 



(25) 



This is also the form of the supermatrix Q used to decouple the interaction via Hubbard-Stratonovich transfor- 
mation; 



Z{u}b,ujf) 



exp 



1, 



-STr 



TV 



Q'+lVN{Q--iu)^a-^a 



(26) 



The integration is performed in the space of matrices Q of the form (25) which is equivalent to the linear 
constraint 

= -Q^- (27) 

More precisely, the integral is taken along the real subspacc in the; complex space (27) where Qb is real, Qf is 
purely imaginary, and X = X* (for convergence of the integral). 
After integrating in we arrive at 



Z{wb,ujf) = DQ [SDet(Q - -^w)] exp U-STr^ j . 



(28) 



Since we are interested in small energy scales (of the order of several level spacings from zero), we can expand 
the action to terms linear in oj and obtain 



Z{ujb,ojf) = J DQ exp ^-ySTr 



r/n2 



^ + lnQ 



+ STr|c2;Q-i 



At large N, the integral is determined by the saddle points of the action 

91 
2 



SoiQ) = STr 



InQ 



By varying the action, the equation of the saddle-point manifold is 

Q' = -l. 

By deforming the integration contour onto the saddle-point manifold, the integral reduces to 



Z{ujb,ujf) = J^DQ exp (^-^So{Q) - STr|a>Q 



(29) 



(30) 



(31) 



(32) 



(the transversal directions do not contribute to the integral because of the supersymmetry) . 

The contour of integration T on the saddle-point manifold should be determined from the condition that the 
original contour of integration can bo deformed onto it without making integrals divergent (sec also discussion 
of this procedure in [1]). For the convergence of the integral (19) over the bosonic components of tp and , 
the energy lob must have an infinitesimal imaginary part Imws < 0. Then, for the convergence of the integral 
(32), the matrix Q must satisfy Im Qb > at infinity on the contour T [Qb is the bosonic diagonal element, 
as shown in (25)]. Besides, the contour F must be compact in the fermionic and non-compact in the bosonic 
sector (see, e.g., [1, 2, 23]). It is shown in [1] that F is the Riemannian symmetric superspace SpO(2|2)/GL(l|l) 
(class CIIL*!!!). 

The key observation, important for taking into account the parity of A^, is that the saddle-point manifold 
F consists of two connected components, which are the images of the two components of the symmetry group 
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SpO(2|2) acting on F. This symmetry group acts on Q by conjugation: Q ^ UQU^^. Explicitly, the two con- 
nected pieces of F may be obtained by rotating by the connected component of the symmetry group SpSO(2|2) 
the two representative matrices 



Qi = 



\ 



( i 



—I 



Q2 = 



\ 



\ 



(33) 



The action (30) is invariant with respect to SpSO(2|2), but acquires an additional shift by 2-Ki between the 
two connected components of the saddle-point manifold. 

It is this property of the supersymmetric action that allows to distinguish between odd and even N in 
the large- A'' limit: in the even-N case (class D, no zero levels), the contributions from the two pieces of the 
saddle-point manifold come with equal signs, and in the odd- A?' case (class B, one zero level) — with opposite 
signs: 

Z{ljb,ujf) = Zi{ljb,i^f) + (-l)^^2(a;B,WF). (34) 
The average spectral density p(a;) may be found as 



1 dZ{ujB,UJF) 

p{uj) = Im ^ 



1 dZ{ujB,U!F) 

-Im ^- 

TT OUlF 



UJB=OJF=UJ — iO 



(35) 



To take the integral over T, we need to parameterize the integration contour: in this section, the param- 
eterization involves the two commuting parameters x,<f> and the two anticommuting ^ and f (we never use 
complex conjugation of anticommuting variables in this paper, and so ^ and ^ should be treated as independent 
variables). The expressions for p{co) in coordinates takes the form: 



p{u;) = Im DQ Qf exp [-Tru){QB - Qf)] ■ 



(36) 



where Qb{x, (f), ^, ^^'^ Qf{x, (j), C) are diagonal matrix elements of Q in a particular parameterization, and 
the measure of integration DQ is 

DQ = ^J{x, 4), Odx d(j) d^. (37) 
27r 

The Jacobian J{x, (j), ^, f ) may be found from expressing the invariant metric STr(rfQ)^ in coordinates 



and taking its superdeterminant 



STr(rfQ)^ = g'^dxidxj 



J({arO) = (SDet5'^y/2. 



(38) 



(39) 



In parameterizing the saddle-point manifold we use the usual trick of splitting the rotation of the supermatrix 
Qi into the two rotations by even and odd generators of the supergroup [2] . Namely, we parameterize 



Q = u^QzU7\ 



(40) 



where Qz is obtained from Qi or Q2 by even rotations parameterized by x and (j) (and without mixing between 
boson-bosonic and fermion-fermionic blocks), and 



C/^ = exp(A), 



(41) 



where A is an odd infinitesimal rotation linear in ^ and ^. 

Supersymmetric calculations of this sort often lead to singularities in superintegrals which need to be resolved 
by properly taking into account boundary terms (see e.g. [2, 1]). In this paper I avoid such singularities by an 
appropriate choice of parameterization of the odd rotation C/^. 
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We shall also employ the symmetry relating the two components of the saddle-point manifold. Namely, 
conjugation by the matrix 



/ 1 
1 



V 



1 

0/ 



(42) 



transforms Qi into Q2 (Q2 — T^^QiT) and the two components of the saddle-point manifold into each other. 

Thus we first parameterize the component generated by Qi, and then the parameterization of the other 
component may be obtained by applying the operator T. 

The even rotations of Qi may be parameterized as 



Qz 



cosh.T — e^'^sinha; 




e~"^sinha; —cosh a; 






-1 




1 



(43) 



(with X e [0; +00), <^ e [0; 2-k]). 

The symmetry (27) of the matrix Q imposes a similar symmetry constraint on the matrix A in (41). The 
latter constraint admits four independent parameters in the boson-fermion and fermion-boson sectors of A. 
However, when acting on Qz, only two of them are independent. At this stage we have a freedom of choosing 
two of the four infinitesimal rotations for our parameterization. The final result does not depend on our choice 
(provided the Jacobian (39) is non-degenerate), but a good choice of parameterization may considerably simplify 
the calculation. 



We choose 



which leads to 





A = 


/ 


e 






e 


} 















°- \ 


V 








i + K? 





1 + 1^"^ / 



(44) 



The Jacobian calculation may be simplified using the simple algebraic identity [2] 

ds' = isTr(dQ)2 = ^STr(dO,)2 + isTr[Q„ ^t/^l^ + STv{5U^ [Q.,dQz\), 



(45) 



where 5U^ = U^^dU^. 

After some calculation, we find for the parameterization chosen 



ds'^ = dx^ + sinh^ x dcjP' - [4 (cosh a; + 1) + 2 sinh^ x i^]d£, d£, + 2i sinh^ x{^d^ + ^ d^)d(f> 



and 



Also, by a direct calculation. 



WBl 

Qfi 



- 1 X 



z [cosh a; — (cosh a; -|- 1)^^], 
i[l + (cosha; + l)^^]. 



(46) 

(47) 
(48) 

(49) 
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Using the operator T to relate the two connected components of the saddle-point manifold, we find for the 
second component 



Qb2 = Qbi, 

Qf2 = -Qfi, (50) 



and the Jacobian obviously remains the same (48). 
As a consistency check, one may verify that 



Z±{lj,lo) = j ^J(a;,^,^,Oc^a;d^d^rf^ (exp[-7rw(Qsi - Qfi)] ±exp[— 7ra;((3si + Qiri)]) = 1 (51) 



(up to a sign). 

Now the calculation of the integral (36) is easily done: 



Pi{u)) = Re / dx i d(j) i d^d^ — tanh - Tl + (cosha; + 1) 
Jo Jo J 47r 2 '- 

xexp(^- i7ra;[(cosha;- 1) - 2(cosha; + 1)^^] ) = ^6{x) + 1, (52) 



f°° f^'^ f - 1 X 

Re dx d(t) I d^d^ — tanh-[l + (coshx + 1)^^] x 
Jo Jo J A-K 2 

/ . r , 1 <-/ ^ sin(27ra;) , „, 

xexp(-j7ra;[cosha; + lj) = -(5(a;) -. (53) 



[all the calculations are performed up to an overall sign]. The ^-function terms are obtained from imaginary 
1/iu! terms by shifting uj to the lower half-plane a; — > w — iO. 

Combining these results with proper signs, we arrive at the final expressions (3) and (4). The asymptotic 
value p{ijj — > cxd) = 1 proves that A given by (16) is indeed the average level spacing. 

Note that pi(w) appeared in Ref. [5] as the spectral density in class BD (which is the average of B and D). 

4 Classes Dlll-even and Dlll-odd 

For classes Dlll-even and Dlll-odd, the calculation is similar to that of the previous section. The saddle- 
point manifold again consists of two connected components, and taking their contributions with different signs 
distinguishes between odd and even matrix dimension. 

The ensembles Dlll-even and Dlll-odd are symplectic (have /? = 4). In the calculation of the average 
spectral density in these ensembles, the matrix Q has dimension 4+4 and belongs to a (8+8)-dimensional linear 
space. The saddle-point manifold is (4-|-4)-dimensional. 

The ensembles Dili are defined as consisting of 2A'' x 2A'' matrices 

H = i(l\ (54) 



—H\^ 

where Hi and H2 are real N x N antisymmetric matrices (Hf = —Hi, iJj = — -^2)- Depending on whether A'' 
is even or odd, this defines the ensemble £>III-even or £)III-odd, respectively. The matrix elements of Hi and 
H2 are assumed to be distributed independently with a Gaussian distribution, and produce the following pair 
correlation function for the matrix elements of H: 

{Hai,bjHa'i',h'3') = V^{Sii'Sjf - {-iy~^^5u'Sjf) {Sab'Sba' - Saa'hb'), (55) 

where the indices i, j take values 1 or 2 and distinguish between the two A-dimensional sectors in the 2^- 
dimensional linear space, and ^ = f ^ ) "1-2" space. 
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We express energy in the units 



A = 



:Trv 



(56) 



(as the resuh of the calculation, this is the average level spacing). 

The space of \l/-vectors needs to be doubled. Instead of a single vector we introduce a pair of vectors 
and ^'2o (here a takes values 1,...,N): 



la 



IplF 
tp2F 

4f 

\4f) 



2a 



/ V'2B \ 
-IplB 

^Ib 

V'2F 
-tjjlF 

^If 
\-4fJ 



/ 4b \ 
-^Ib 

1p2B 

4f 
-■^If 

fplF 
-■ 02F, 



2a 



/ 4b ■ 

4b 

~tp2B 
-fplB 

V'2F 



(57) 



The two sets of vectors 4*, "if are necessary to reproduce the four terms in the interaction induced by (55). 
The corresponding supermatrix Q has the form 



Q 



( Qb 


Xb 


-Yb 





0-1 


-02 


Pl 


-P2 




Xb 


-Qb 





Yb 


-a2 


-U\ 


-P2 


-pl 




Yb 





-Qb 


Xb 


pl 


-P2 


o\ 


-(J2 







-Yb 


Xb 


Qb 


-P2 


-Pl 


-a 2 


-CTl 




<y\ 


0-2 


-Pi 


-P2 


-Qf 


-Xf 





-Yf 




0-2 


-(Tl 


-P2 


Pl 


-Xf 


Qf 


-Yf 







Pi 


P2 


-ai 


-0-2 





-Yf 


Qf 


-Xf 




\ P2 


-Pl 


-CT2 


CTl 


-Yf 





-Xf 


-Qf 


) 



(58) 



FB,PH,12 



Equivalently, Q may be described as obeying the two linear constraints: 



7l2<37l2^ = -Q, 



712 = 



-1 



(59) 



and 



IphQIph 



Q' 



JPH = 7 



1 






-1 











1 




I—' 






FB,PH 



(60) 



where FB and PH indices specifies that the operator acts in the Fermi-Bose and "particle-hole" spaces (the 
doubling of dimension by combining ip and ?/'''' in a single vector '3/), and "1-2" denotes the space corresponding 
to the two iV-dimensional sectors in the original Hamiltonian (54). 

Similarly to the previous section, we find for the generating function Z{u)B,tJF) 



Z{ujb,u}f) 



I 



exp I -^STr 



-j^^ia-^iaCO + -(*ia*ia) 



= j DQ J D'i exp I -^STr 



y<3'+*YY(Q-;^'^)*ia*i 



DQ |^SDet(g - — w) 



-N/2 



exp I -ySTr— 



(61) 
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where 



w = 



U>B 






UJf 



1 

-1 



PH 



1 

-1 



(62) 
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At small energies lvb, <^f this leads to the formulas (29)-(32), albeit with the new definitions of Q and lo. 
The saddle-point manifold consists of the two connected pieces represented by 



Qi = 



—I 



—I 



\ 



\ 



—I 



( i 



Q2 



—I 



—I 



\ 



\ 



—I 



—I 



—I 



—I 



(63) 



Similarly to the procedure described in the previous section, first the supermatrices Qi, Qi are rotated by even 
symmetry-group generators. These rotations do not mix bosonic and fermionic components, i.e. the matrix Qz 
contains only boson-boson and fermion-fermion blocks: 



(64) 



We shall use the following parameterization of these blocks: 

/ icosh^B nisinh^B (n2 — ins) sinh^B 



i cosh 6b 

rii sinh 9b 
(n2 + ins) sinh^B 




ni sinh^B 
— icosh 9b 


-(n2 -h ins) sinh ^B 





—i cosh^B 
ni sinh^B 



-(n2 — in^) sinh 6'b 
ni sinh 9b 
i cosh 9b / 



(65) 



PH,12 



where (ni, n2, ns) is a vector of a real two-dimensional unit sphere (n{ -|- n| -|- 77.3 = 1). The boson-boson block 
is the same for the two sectors of the saddle-point manifold. 

The fermion-fermion blocks for the two components of the saddle-point manifold are: 












-ie 




































-i COS 9p —isui9p 

-ism9F icos9p 

icosOp —ism9p 

— isin^iT —icos9p, 



(66) 

The parameter range is 9b S [0; -l-oo), 9p G [0; 27r], and the vector n runs over the two-dimensional unit sphere 
52. 

Like in the previous section, we first do the calculation in the first component of the saddle-point mani- 
fold, and then obtain the answers for the second component by using the symmetry operator T mapping one 
component (generated by Qi) onto the other (generated by Q2)- One possible choice of such a matrix T is 



/ 1 



1/2 -i/2 -i/2 1/2 

-i/2 -1/2 1/2 i/2 

-i/2 1/2 -1/2 i/2 

1/2 i/2 i/2 1/2 ) 



(67) 



FB,PHA2 
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Returning to the parameterization for the first component of the saddle-point manifold, the matrix A involved 
in the odd rotation (40,41) is chosen as follows: 

/ 



A = 













V 





\ 










— V 


















V 


I 
















-I 


V 






















































—1/ 














v 













(68) 



FB,PHA1 



This matrix satisfies the flatness condition [A, rfA] = 0, and this leads to 6U^ = ^dU^ = dA. Using the 
algebraic identity (46), together with (38), (39), and 



DQ = 



:J{9b, Of, n, ^, v) dOs dOp d^nd^ dS^du dv, 



one finds after some calculation the explicit form for the invariant measure in the coordinates chosen: 



DQ 



1 



167r2 



-2i0F 



sinh 6b d^B dOp d nd^d^dudv 



(here d^n is the integration over the solid angle on the unit sphere). 

The explicit expressions for the diagonal entries of the Q-matrix Qb and Qp are found to be 

Qbi = i[cosh0B + e*«-(e^-Pi.)], 
Qp, = -ie'^^{^^ + i?iy). 

Using the operator T defined in (67), for the second component of the saddle-point manifold we find 

Qb2 = Qbi, 

Qf2 = i [cos Op + 2^^j>j/e'^^ + cosh enivv- 

+ini sinh0_B(i^C ~ ^0 ~ ("2 + in^)^^, sinh 0b + (n2 — ma)!/^ sinh 0b]- 

After some calculation, one verifies the normalization: 

Zi{(jj, ±w) = j dQ exp[-27ra;((3Bi T Qfi)] = 0, 

Z2{lj, ±uj) = j dQ exp[-27rw(QB2 T Qf2)] = 1- 
The density of states is found in terms of Bessel functions: 

f I IT 

pi{uj) = Im / dQQpieyi^[-2-KU){QBi - Qfi)] = -S{uj) - -Ji{2'kuj) 
in the first sector, and 

P2{ijj) = dQQB2exp[-27ra;(QB2 - Qf2)] = ^5{ijj)+-k^w [J((27ra;) Jo(27ra;) + Ji^(27ra;)] , 



(69) 



(70) 



(71) 



(72) 



(73) 



(74) 



(75) 



where in (74) and (75) we assumed w > and are careless about the overall sign of the answers. 

Taking these contributions with proper signs, we obtain the answers (5) and (6). Remarkably, in (6), the 
contributions of Pi{lo) and of P2{'-^) cancel each other to the third order at small uj, producing the correct 
behaviour of the total density of states p{io) oc lo^ . 
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5 Class yllll (chiral unitary) 



The supersymmetric calculations for the three chiral random-matrix ensembles differ from those for classes B-D 
and Dili in that the broken symmetry of the saddle-point manifold is not the discrete Z2, but the continuous 
GL(1). The representations of GL(1) are enumerated by the integer winding number whose absolute value 
equals the number of zero levels in the random-matrix ensemble. 

The chiral unitary ensemble considered in this section has (3 = 2 (unitary bulk statistics). In the calculation 
of the average density of states, the supermatrix Q has the block form (in the "1-2" space) 



Qi 
Q2 



(76) 



where Qi and Q2 arc (l+l)-dimcnsional supermatrices without linear constraints. Thus the linear space of 
matrices Q is 4+4-dimensional, and the saddle-point manifold has dimension 2+2. 

The (Gaussian) chiral unitary ensemble (class ^4111 in Cartan notation) consists of the matrices of the form 



H 



H 




(77) 



where if is a rectangular matrix pxq with complex matrix elements. The matrix elements of H have independent 
Gaussian distributions: 



dP{H) oc Yl 

exp I — - 

a,b V 



dReHab dlmHab- 



(78) 



The spectrum of such a matrix H consist of N = mm{p, q) pairs of eigenvalues ±Ei and of m = |p — zero 

eigenvalues. In this paper wc arc interested in the average density of states near zero in the limit of large 
matrices N 00 while keeping the number of zero levels m fixed. The average level spacing near zero in the 
large- A?^ limit is 



A = 



2VN 



(79) 



Following the standard procedure, we introduce p-component superfields and ^'i, and g-component 
superfields ^'2 and ^'2: 



*26 



i>2F 



*26 



i'lp/b 



a= l,...,jj, 
h=l,...,q. 



(80) 
(81) 



The generating function (17) takes the form: 



Z{ujb-,u>f) 



L>* exp -STr 



2^/N 



(*la*la + *26*2b) OJ + (*la*la)(*26*26) 



= j D{Qi, Q2) J exp (-STr [NQ1Q2 + i^/N (Qi - ^2b^2b + i^N (Q2 - ^w) ^-la^-ia] ) = 

(82) 



/ 



i3(0i,Q2) [SDet(gi-— w) 



where the original integration contour in Q is at Q2 =Q\-, and 



SDet(g2 - ^w)] " exp {-N STr Q1Q2) 
-nt 



At small energies, expanding the supersymmetric action to terms linear in u)b,f, we find 



Z{u)b,ijJf) 



DQ [SDetQi]"exp [-NSo{QuQ2) + -S'IVw(Qr' + Q^') 



(83) 



(84) 
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with 



5o((3i, (32) = STr (Q1Q2 + In Qi + In Qs) . 



(85) 



If Qi and Q2 are combined into a single supermatrix Q according to (76), this action coincides with the standard 
form (30). The saddle-point manifold T is given by the condition 



Q1Q2 = -1 



(86) 



[which is equivalent to = —1], and the generating function may be written as the integral over the saddle- 
point manifold 

. . .1 (3,^ 



Zm{LOB,uJF) = J [SDetQiJ^'exp ^--ST^w(gi + O2) 



A paranictc;rization and the calculation of the integral was previously performed in Ref. [6], and I outline 
their calculation here for completeness. The matrices Qi and Q2 are parameterized as 



Qi = —Q2 ^ 

















exp 










, 



(88) 



The invariant measure deduced from the metric ds"^ = STr{QiQ2) leads to the trivial Jacobian J = 1 and to 
the integration measure 

DQ = ^dxd<pd£,dl (89) 

From the parameterization (88), the diagonal elements of the matrices Qi and Q2 are easily computed. The 
generating function (87) involves the average of the diagonal elements of Q\ and Q2: 



Qbb = 2 [(<3i)bs + (Q2)bb] = i cosha; 



Q 



FF 



[{Qi)ff + {Q2)ff] = icosc 



Also, the m-dependent prefactor in (87) is a plane wave generated by 

SDetQi = e^-'"^. 

The normalization can be verified by computing the integral 

Zm{oJ,oj)= [ ^dx d(t) d^ di e™^^-''^) 
J 27r 

The average spectral density is calculated from (35) as the integral (again, up to an overall sign): 

p„(a;)=Re / ^ dx#d^d^e'"(^-''^' cosh 
which after some algebra produces the result (7). 

















exp 


—iiTUJ ^coshx 




+ coscj) 




)] 


= ±1 



(90) 



(91) 



(92) 







—inoj ^cosha; 










Hi 


exp 


Hi 


+ COS(/) 


Hi 


)1 



(93) 



6 Class BDI (chiral orthogonal) 



Similarly to the three Wigner Dyson random-matrix ensembles, the chiral ensembles form the three classes: 
unitary, orthogonal, and symplectic, depending on the structure of the matrix H in the block form (77). In 
the chiral orthogonal ensemble (class BDI in Cartan notation), the matrix H is real. In this ensemble, the 
bulk repulsion of the levels corresponds to the orthogonal regime: [3 = 1. The supermatrix Q involved in the 
calculation of the average spectral density has the block form (76), but now the matrices Qi and Q2 have 
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dimensions 2+2 each, with one linear constraint, so the dimension of the linear space of all the matrices Q is 
8+8. The saddle-point equation selects the saddle-point manifold of dimension 4+4. 

The entries of the matrix H in (77) are assumed to be real, with independent Gaussian distributions: 



dP{H) oc Y[ exp 



a, 6 



\Hab\- 

2^2 



dHab- 



(94) 



As in the chiral unitary ensemble, the spectrum consists of N pairs of opposite energies, and of m = |p — 
zero-energy levels. The average level spacing is given by the; same cixprcssion (79) [note however a difference 
between the definitions of v in (78) and in (94), depending on whether Hab are complex or real]. 

To account for the matrix elements Hab being real, we need to double the dimensions of the superfields 
and ^j. The p-component superfields ^'i and are defined as 









(4b\ 








^IB 

4f 




a 





a = 1, . . 



(95) 



and similarly, the g-component superfields '^2 and '^2'- 









( 4b \ 








i'2B 




b 


V -i^2F / 



1,.. 



(96) 



Repeating the steps of the derivation (82), we arrive at the following expression for the generating function [in 
place of (84)]: 



Z^{ojb,ujf)^ y"DQ [SDetQi]"/'exp So{Qi, Q2) + ^STr 6j{Q^' + Q2') 



(97) 



where uj and So{Qi,Q2) are given by the old expressions (83) and (85). However, the matrices Qi and Q2 are 
now two times bigger. Each of them has the explicit form: 



Qbi 








X, 


Qbz 


P^ 






pi 


Qfi 





Pi 


(Ti 





Qfi , 



Qi 



Equivalently, this form of the matrices Qi may be described by the linear constraints 

7Qi7-' = Ql 

where 



(98) 



/O 1 
1 



1 

V -1 oy 

The saddle-point manifold T is determined by the condition (86), and the generating function is expressed as 



(99) 



(100) 



Zm{LOB,tOF) 



DQ [SDetQi]'"/^ 



exp 



--STra;(Qi + Q2) 



(101) 
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Ill the present section we choose a shghtly different form of parameterization than in the previous one. 
Namely, parameterize 

Qi = UiQ,iU^\ Q2 = U2Q,2U^\ (102) 

where the matrices Qzi and Qz2 contain only boson-boson and fermion-fermion blocks (contain only even 
rotations), and the matrices Ui and U2 contain only odd rotations. The explicit form of this parameterization 
is as follows: 



Qzi = i 



1 e^cosh^ e^+'2'sinhi9 






( 


e ^cosh^ — e ^+*^sinh^ 






, Qzi = i 




_ g-X-SJ/ gjjjjj g-X ^.Qgjj ^ 






e*"^ 






e-"^ 


V 


e'"^ ) 




\ 




e-''^ 



Ai = 



C/i = cxp(Ai 
-A i 



\ 



U2 = exp(^2), 

A2 = 



i A 



After some calculation, the Jacobian is found to be 

J = sinh6»e2(^-^*), 

and therefore the measure of integration is 

DQ = sinh e e^^-"P) dx dy dO dcp d^ d^ d\ dX. 

(27r)2 

Finally, from the explicit calculation of the diagonal elements of Qi and Q2- 



(103) 
(104) 

(105) 

(106) 
(107) 



Qbb = ^ [{Qi)bb + {Q2)bb] = i 
Qff = ^ 1{Qi)ff + {Q2)ff] = i 



coshcc cosh 61 - ^e"t'{^^ - XX) 



cos(t)+ -e-^ (cosh 61 [^^ - AA] - sinh6l [e'^^A + e-'^'A^]) 



Now, after verifying the normalization 

ZmiiJ,±uj) = J £»Qe™(^-'^)exp[-7ra;(QsB ±(3ff)] = 1, 
we compute the average density of states as the following integral: 

Pm{uj) = ReJ £»ge™(^-**) Qbb exp[-7ra;(gsB + Qff)], 



. (108) 



(109) 



(110) 



[where DQ, Qbb, and Qff are defined in (107) and (108)]. After some algebra and manipulations with Bessel 
functions, this produces the result (8). 



7 Class CII (chiral symplectic) 



The last symmetry class considered in this paper is CII in Cartan notation: the chiral symplectic one. It consists 
of the matrices H of the block form (77), where the matrix H has an internal 2x2 structure: 



H = 



a b 
a' 



(111) 
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where a and b are rectangular p x gr matrices [the matrix H tlnis has the dimensions 2p x 2q, and tlic Hamihonian 
(77) has the dimension 2{p + q)]. The spectrum of such a Hamiltonian consists of the N = min(p, q) pairs of 
doubly degenerate levels at opposite energies ±Ei, and of 2m = 2\p~ q\ zero-energy levels. Since each level has 
degeneracy two (or a multiple of two), we divide the density of states by two for the purpose of level counting, 
and count each degenerate level as a single one. 

As in the two previous sections, the supermatrix Q is of the form (76). Similarly to the chiral orthogonal 
case, the matrices Qi and Q2 in this calculation have dimension 2+2, with one linear constraint. The dimension 
of the linear space of matrices Q is thus 8+8, and the dimension of the saddle-point manifold is 4+4. 

The entries of the matrix (111) are assumed to be normalized as 



dP{A) (X Yl 

exp 

a,b 



\aab\ 



dRettab dlT^ciab JJoxp 



a.b 



\bab 



dRebab dlrabab- 



With this normalization, the average level spacing is 

A = 



'2N 



(112) 



(113) 



Similarly to the fill ensemble discussed before, we need to additionally double the set of the superfields to 
account for the symplectic matrix structure. Namely, in addition to the "1-2" sector, distinguishing between 

the sizc-p and sizc-g blocks in (77), we introduce the "particle- hole" (PH) sector referring to the two sectors 
inside the matrix A (111). Thus we arrive to the four pairs of superfields \l/pi, \E'pi, ^ P2, *P2j 



Pla 



Pla 



/ IpPlB 

--(pPlF 



Hla 



( <IB \ 

-IpPlB 

\-^pif) ^ 



Hla 




and similarly for the sector-2 fields ^'p2, *P2, *//2, ^ H2- 
This results in the following form of matrices Qi: 



Qi 



a=l, 

(114) 



/ Qbz 





-CTi 


Pi 


\ 





Qsi 


Pi 








-Pi 


Qpi 


X^ 




\ Pi 


CTi 




Qpi 


J 



(115) 



FB,PH 



which differs from the chiral orthogonal ensemble by interchanging the bosonic and fermionic sectors (this 
duality was already described in Ref. [1]). 

Performing the standard steps of the derivation [similarly to (82)] , we arrive at the following answer in the 
saddle-point approximation: 



Zm{0JB,0JF) 



DQ [SDetQi]™exp [--STrw(Qi + Q2) 



(116) 



[note the differences from the orthogonal result (101)!] 

Finally, the calculation may be performed using the parameterization obtained from that of the previous 

section by interchanging bosonic and fermionic sectors. Note that under this interchange, the compact variables 
become non-compact and vice versa. Explicitly, eqs. (103), (105)-(107) are replaces by 



Q 



zl 



















e*"* cos 


e^(0+y) sin 6* 




-e^i't'-y) sin6» 


e**cos6l , 



Q 



z1 



e-^ 






\ 


e"^ 










e-^-^cos^ 


-e.-^{<l>-y) sin 61 






(,-i{<l>+y) sine 


e"'"^ cos 6* 


/ 



(117) 
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V 



-A 



e A 
VA C 



DQ = sin e e2(^-''^) dx dy d0 dcj) d^ d^ dX dX. 

(27r)2 



(118) 
(119) 



Eq. (108) is replaced by 



Qbb = 2 [(<9i)bb + {Q2)bb] = i 
Qff = \ [{Qi)ff + {Q2)ff] = i 



coshx + -e"^ (cos6» [^^ + AA] + sin6' [e*^A^ + e'^^A^]) 



cos0cos^--e ''(^^ + AA) 



Taking the integral 



=ReJ 



DQe 



2m{x—i(j)) 



Qbb exp[-27rw((5BB + Qff)], 



(120) 



(121) 



we obtain the final result (9). [Note that in (121) we divided the density of states by two to prevent double 
counting of the doubly degenerate states.] 



8 Zero levels and reduced super symmetry of the action 

After presenting the calculations of the average density of states in the five random-matrix ensembles with zero 
levels, in this section we discuss the specifics of the supersymmetric method due to the zero levels. I do not 
present here a consistent mathematical analysis of this problem, leaving it for future study. Instead I only 
summarize the common features of the above calculations specific for the ensembles with zero levels. 

The standard supersymmetric procedure to calculate spectral correlation functions in any random-matrix 
ensemble starts with introducing bosonic and fermionic fields i])b and i]jp [1, 2]. Integrating over the Gaus- 
sian disorder produces a four-term interaction. This interaction is then decoupled via Hubbard-Stratonovich 
transformation by a supermatrix Q whose dimension is independent on the matrix size A'' in the original 
random-matrix ensemble. Integrating over the superfields ('0b, V-'f), one arrives at an effective action for the 
supermatrix Q. The supermatrix Q obeys certain linear symmetry relations and thus belongs to a linear su- 
perspace L depending on the symmetries of the original random-matrix ensemble. In the superspace L, there 
acts a supergroup G inherited from the supersymmetry mixing the bosonic and fermionic fields V'b and V'f. 
Namely, an element G G is a supermatrix acting on Q by conjugation: 



g:Q^ gQg~ 



(122) 



The group G then depends on the symmetries of the matrix Q and hence on the symmetry of the random- matrix 
ensemble. 

The effective action S{Q) may, at small energies w, be expanded to the linear in w term: 



S{Q) = 5„=o(Q) +wSTrAQ-i, 



(123) 



where A G L is a particular supermatrix with A^ = 1. 

In ensembles without zero eigenvalues, the zero-energy action S^=o{Q) is invariant with respect to the 

supergroup G. The zero-energy action St^=o{N) scales linearly with TV, and in the large- A'^ limit the integral over 
Q is determined by the saddle points of S'aj=o {Q) ■ The saddle-point equation, with the appropriate normalization, 
reads 

= -1. (124) 

This equation is solved by Q = «A, as well as by any matrix obtained from iA by G- rotations (122). The matrix 
A is invariant under a subgroup H C G, and the saddle-point manifold F is the quotient G/H. 
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Table 2: Supergroups and superspaces involved in the spectral density calculations 



Class 


L 


G 


H 


r = G/H 


Go 


G/Go 


B-D 


osp(2|2) 


SpO(2|2) 


GL(1|1) 


SpO(2|2)/GL(l|l) 


SpSO(2|2) 


Za 


Dili 


osp(4|4)/osp(2|2)eosp(2|2) 


SpO(2|2)xSpO(2|2) 


SpO(2|2)* 


SpO{2|2) 


S[SpO{2|2)xSpO{2l2)]** 


Z2 


AIII 


gi(i|i)egi(i|i) 


GL(1 l)xGL(l|l) 


GL(1|1)* 


GL(1|1) 


S[GL(1 l)xGL(l|l)]'* 


GL(1) 


BDl 


[gl(2|2)/osp(2|2)]e[gl(2|2)/osp(2|2)] 


GL(2|2) 


OSp(2|2) 


GL(2|2)/OSp(2|2) 


SL(2|2)*** 


GL(1) 


CU 


[gl(2|2)/osp(2|2)]e[gl(2|2)/osp(2|2)] 


GL(2 2) 


SpO(2 2) 


GL(2l2)/SpO(2|2) 


SL2(2|2)*** 


GL(1) 



* H is diagonal in G: H = {{g,g)} in classes Dili and Alll 

** S[H X H] denotes here the subgroup {(gi, g2)|SDet(/i = SDetga} 

*** SL(2|2) = {ge GL(2|2)|SDetg = 1}; SL2(2|2) = {g 6 GL(2|2)|SDetg = ±1} 



More precisely, F is a Rieniannian (real) supermanifold in G/H, which makes it a Riemannian symmetric 
superspacc as defined in Ref. [1]. This real submanifold should be determined geometrically from deforming 
the real integration subspace in L while keeping the integral convergent. A good geometric understanding of 
this contour deformation still needs to be developed, but the rule of thumb for choosing the real integration 
manifold in F is to take the bosonic sector non-compact and the fcrmionic one compact [this choice also provides 
a metric of a definite sign on F]. 

For the random-matrix ensembles with zero eigenvalues, the invariance properties of the effective action 
Suj=o{Q) are modified. In this case, Suj=o{Q) is invariant with respect not to the whole supergroup G, but only 
with respect to its normal subgroup Go- The subgroup Go must contain H, and the factorgroup G/Gq is an 
abelian group (an ordinary group, not a supergroup). The exponent exp[S'^=o('5)] transforms according to one 
of its (one-dimensional) representations. The degree of this representation equals the number of zero eigenvalues 
in the random-matrix ensemble. 

In the ensembles B-D and Dill, the group G/Gq is discrete Z2, and the two representations of Z2 correspond 
to the ensembles with and without zero eigenvalues (odd A'' and even N, respectively). Specifically, the action 
Soj=o{Q) has the form 

S^^oiQ)=NSoiQ), (125) 

where So{Q) gets incremented by in under the action of the generator of G/Gq = Z2. Hence, exp[S'a,=o(Q)] 
transforms according to the even/odd representation of Z2 for even/odd N. 

In the chiral random-matrix ensembles ^III, BDl, and GII, the group G/Gq is the continuous GL(1), with 
its representations labeled by the integer "winding number" m. The absolute value of m equals the number of 
zero eigenvalues in the random-matrix ensemble. The action Sui=o{Q) in the chiral ensembles [which includes 
the logarithm of the pre-exponent in (87), (101), (116)] is of the form 

Su.=o{Q) = NSoiQ) + mSiiQ), (126) 

where So{Q) is invariant under G, and Si{Q) produces phase shifts under G/Gq = GL(1). The exponent 
exp[S't^=o(Q)] then transforms as the representation of GL(1) of degree m. 

The summary of the "building blocks" of the supcrsymmetric calculations in this paper (the calculation of 
the average spectral density) is presented in Table 2. This table is compiled using the results of Ref. [1] and the 
calculations in the previous sections. The definitions of the supergroups involved in this table may be found 
in Section 2. To keep track of the bosonic (non-compact) and fcrmionic (compact) sectors of the supergroups, 
we use the notation OSp for the orthosymplectic supergroup with the orthogonal part in the bosonic, and the 
symplectic part in the fermionic sectors. In the opposite case of orthogonal fermionic and symplectic bosonic 
sector, we denote the same supergroup SpO. It is important for reducing the action symmetry in ensembles 
B-D and Dili, that SpO(2n|2n) has two disconnected components with superdeterminants 1 and —1 [the 
former of them denoted as SpSO(2n|2n)]. Incidentally, in the ensembles G and GI, dual to B-D and Dili by 
interchanging fermionic and bosonic sectors, in the supergroup OSp(2n|2n), with the non-compact orthogonal 
sector, the second component (with SDet = — 1) plays no role and should be disregarded as it always corresponds 
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to a divergent integral. The supergroup GL(n|n) is not simple cither. Firstly, it has a one-dimensional center 
consisting of scalar matrices. Secondly, it has a normal subgroup SL(n|n) consisting of matrices with unit 
superdeterminant. This latter reduction of the STipergroup GL(n|n) is crucial for the symmetry classification in 
the case of the chiral ensembles (the last three lines in Table 2). 

Finally, it is worth mentioning that the conclusion about the reduced supersymmetry of the zero-energy 
effective action 5*^^=0 (Q) may be extended to higher-order correlation functions involving averaging several 
Green's functions [the average spectral density requires averaging only one Green's function]. In Ref. [1] a 
general procedure of calculating correlation functions of arbitrary order (with the number of zero levels m = 
0) was described, and the saddle-point manifold (a Riemannian symmetric superspace) F was found to be 
always reducible for ensembles admitting zero levels. Thus the extension of the calculation of the present 
paper to higher-order correlations is straightforward (however, explicit parameterization and integral evaluation 
immediately becomes much more complicated). 

I thank N. Nekrasov, P. Ostrovsky, and M. Skvortsov for teaching me different aspects of supersymmetry, 
and Swiss National Foundation for financial support. 

References 

[1] M. R. Zirnbauer, "Riemannian symmetric superspaces and their origin in random-matrix theory", 
J. Math. Phys. 37 (1996), 4986 [math-ph/9808012]. 

[2] K. Efetov, "Supersymmetry in disorder and chaos", Cambridge University Press, Cambridge 1997. 

[3] S. Helgason, "Differential geometry. Lie groups, and symmetric spaces", Academic Press, New York 1978. 

[4] M. Caselle, "A new classification scheme for random matrix theories", cond-mat/9610017. 

[5] M. Bocquet, D. Serban, and M. R. Zirnbauer, "Disordered 2d quasiparticles in class D: Dirac fermions with 
random mass, and dirty superconductors", Nucl. Phys. B 578 (2000), 628 [cond-mat/9910480]. 

[6] P. H. Damgaard, J. C. Osborn, D. Toublan, and J. J. M. Verbaarschot, "The microscopic spectral density 
of the QCD Dirac operator", Nucl. Phys. B 547 (1999), 305. 

[7] D. A. Ivanov, "Random-matrix ensembles in p-wave vortices", to appear in "Vortices in unconventional su- 
perconductors and superfluids — microscopic structure and dynamics" , Springer, 2001 [cond-mat /0103089] . 

[8] J. Verbaarschot, "Spectrum of the QCD Dirac operator and chiral random matrix theory" , Phys. Rev. Lett. 
72 (1994), 2531. 

[9] J. J. M. Verbaarschot and T. Wettig, "Random matrix theory and chiral symmetry in QCD", hep- 
ph/0003017. 

[10] M. L. Mehta, "Random Matrices", Academic Press, Boston 1991. 

[11] T. Nagao and K. Slevin, "Laguerre ensembles of random matrices: Nonuniversal correlation functions", 
J. Math. Phys. 34 (1993), 2317. 

[12] A. Altland and M. R. Zirnbauer, "Nonstandard symmetry classes in mesoscopic normal-superconducting 
hybrid structures", Phys. Rev. B 55 (1997), 1142 [cond-mat/9602137]. 

[13] J. J. M. Verbaarschot and I. Zahed, "Spectral density of the QCD Dirac operator near zero virtuality", 
Phys. Rev. Lett. 70 (1993), 3852. 

[14] T. Nagao and K. Slevin, "Nonuniversal correlations for random matrix ensembles", J. Math. Phys. 34 
(1993), 2075. 

[15] P. J. Forrester, "The spectrum edge of random matrix ensembles", Nucl. Phys. B 402 (1993), 709. 



21 



[16] A. V. Andreev, B. D. Simons, and N. Taniguchi, "Supersymmetry applied to the spectrum edge of random 
matrix ensembles", Nucl. Phys. B 432 (1994), 487. 

[17] J. Verbaarschot, "The spectrum of the Dirac operator near zero virtuahty for Nc = 2 and chiral random 
matrix theory", Nucl. Phys. B 426 (1994), 559. 

[18] T. Nagao and P. J. Forrester, "Asymptotic correlations at the spectrum edge of random matrices", 
Nucl. Phys. B 435 (1995), 401. 

[19] T. Wettig, "Random-matrix theory in quantum chromodynamics and on the lattice", Habilitationsschrift, 
Ruprecht-Kar Is- University Heidelberg, April 1998. 

[20] B. Klein and J. J. M. Verbaarschot, "Spectral universality of real random matrix ensembles", Nucl. Phys. 
B 588 (2000), 483. 

[21] J.-Z. Ma, T. Guhr, and T. Wettig, "Statistical properties at the spectrum edge of the QCD Dirac operator" , 
Eur. Phys. J. A2 (1998), 87, erratum ibid. A2 (1998), 425 [hep-lat/9712026]. 

[22] B. DeWitt, "Supermanifolds" , Cambridge University Press, Cambridge, 1992 (second edition). 

[23] J. J. M. Verbaarschot, H. A. Weidenmiiller, and M. R. Zirnbauer, "Grassmann integration in stochastic 
quantum physics: the case of compound-nucleus scattering", Phys. Rep. 6 (1985), 367. 



22 



